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3l.  \introduction 

The  purpose  of  this  paper  is  to  study  the  traction  problem  in 
three  dimensional  nonlinear  elasticity  using  geometric  techniques  and 
singularity  theory.  The  first  two  napers  in  the  series  deal  with  dead 
loads  and  configurations  that  are  nearly  stress  free.  As  was  shown  by 
Signorini  "1930]  and  Stoppelli  [1958],  this  problem  has  nontrivial 
solutions.  However,  their  analysis  is  incomplete  for  three  reasons. 
First,  their  load  was  varied  only  by  a  scalar  factor;  in  a  full 
neighborhood  in  load  space  of  a  load  which  has  an  axis  of  equilibrium 
there  are  additional  solutions  missed  by  their  analysis.  Second, 
their  analysis  is  only  local  in  the  rotation  group,  so  additional 
nearly  stress  free  solutions  are  missed  by  restricting  to  rotations 
near  the  identity.  Third,  some  classes  of  loads  with  a  degenerate  axis 
of  equilibrium  were  not  considered.  This  paper  completes  their  analysis 
by  treating  these  questions  as  well  as  stability.  The  complexity  of 
the  answer  is  indicated  by  the  fact  that  near  certain  tv-pes  of  loads, 
we  find  up  to  40  distinct  solutions  that  are  nearly  stress  free.  Our 
constitutive  hypotheses  on  the  stress  tensor  are  'generic';  for  a 
degenerate  stress  tensor  there  can  be  even  more  solutions. 

The  literature  on  this  problem  is  very  extensive,  going  back 
to  Signorini  in  the  1930' s.  Our  primary  sources  have  been  Stoppelli 
[1958],  Grioli  [1963],  Truesdell  and  Noll  [1965],  Van  Buren  [19b8], 

Wang  and  Truesdell  [I9T3],  and  Capriz  and  Podio  Guidugli  [1974], 

However,  none  of  the  references  beyond  Stoppelli  [1958]  proves  any 
of  the  theorems  dealing  with  nontrivial  cases;  i.e.  loads  with  axes 
of  equilibrium.  However,  Grioli  [1962]  is  a  convenient  reference 


for  the  statements. 


The  outline  of  this  first  part  of  the  paper  is  as  follows.  Our 
notation  for  nonlinear  elasticity  and  the  problem  near  a  natural  state 
is  formulated  in  Section  2.  In  Section  3  the  basic  properties  of  the 
astatic  load  are  reviewed  and  developed.  The  problem  is  reformulated 
with  special  reference  to  the  global  aspects  of  the  rotation  group  in 
Section  4  and  introduces  the  bifurcation  equation  which  plays  a  crucial 
role  throughout  the  paper.  Section  5  treats  loads  with  no  axis  of 
equilibrium;  there  are  three  new  features  in  this  section.  First  the 
proof  of  Stoppelli's  results  is  considerably  simplified.  Second,  the 
results  are  global  relative  to  the  rotation  group.  Finally,  the  stabi¬ 
lity  of  the  solutions  is  determined.  The  number  of  solutions  will  be 
classified  by  load  types;  this  classification  scheme  is  explained  in 
Section  6.  (Some  work  related  to  the  "type  classification"  was  given  by 
Ogden  [1977]).  In  Section  7  a  second  order  bifurcation  equation  is  shown 
to  be  a  gradient.  This  consequence  of  Betti  reciprocity  is  basic  to  our 
analysis.  Section  8  gives  a  complete  bifurcation  analysis  of  type  1 
loads  (the  case  considered  by  Stoppelli),  including  a  stability  analysis. 
.Vew  local  and  global  solutions  are  found.  The  final  section  makes 
explicit  the  comparison  with  Stoppelli's  theorem. 

The  second  paper  of  this  series  will  analyse  the  remaining  types 
2,  3  and  4,  using  a  reformulation  of  our  gradient  results,  discuss 
linearisation  stability,  parallel  loads  and  will  give  additional 

connections  with  the  literature. 
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Statement  of  the  Problem 


5  * 

Let  B  C  ]R  be  an  open  bounded  set  with  smooth  boundary  and 

assume  for  convenience  that  0  S  8.  Let  1  <  p  <  s  >  f3/p)  *  1 

and  let  C  be  the  space  of  maps  .-p :  8  !R  that  are  of  class  W 

1  SO 

(so  they  are  C  )  such  that  p(0)  =0,  j  is  a  iV  ’^-diffeomor y'  sm 

onto  its  image,  and  J('v)  >  0,  where  J(j)  is  the  Jacobian  of  . 

For  example,  if  8  ^  is  W^’^-close  to  the  identity  and 

‘;;(0)  =  0,  then  'it  S  C  .  If  Q  is  a  linear  isomorphism  of  IR'’  to 

with  det  Q  >  0,  then  Qop  €  C  as  well. 

Let  points  in  8  be  denoted  X  S  8  and  points  in  TR'^  be 

denoted  x.  Sometimes  we  write  x  =  c(X).  Let  T^8  be  the  tangent 

space  to  8  at  X,  regarded  as  vectors  in  IR'^  based  at  X,  We  do 

not  identify  T^S  and  IR'^  for  conceptual  clarity.  For  c  -  C  ,  let 

F(X)  S  L(T^B,  IR'^)  be  the  derivative  of  i  at  X;  by  standard  abuse 

of  notation  we  write  F(X)  =  D(i(X)  or  74  (X)  interchangeably  and 


L(T,^8,  IR"^)  denotes  the  set 

of  all 

linear  maps 

of 

T^s 

to  IR"^ . 

We  let  FCX)’’'  €  L(1R'',  T^81 

denote 

the  adjoint 

of 

F(X) 

relative  to 

the  Euclidean  inner  product.  Observe  that  F(X]  £  L^(T^B,  IR'^1  ,  the 
linear  transformations  with  positive  determinant,  since 
det  F(X)  =  J(p)(X)  >  0.  We  let  C  =  F^F  (that  is, 

CfXl  =  F(X1^F(X)  €  L(T^8,  7^8)  )  denote  the  Cauchy-Green  tensor. 


We  believe  that  our  results  also  hold  when  6  has  piecewise  smooth 
boundary.  This  program  depends  on  elliptic  regularity  for  such  regions. 
Except  in  special  cases,  this  theory  is  non-existent  and  seems  to 
depend  on  a  modification  of  the  usual  Sobolev  spaces  near  corners. 
However,  for  simple  shapes  like  cubes,  the  necessary  regularity  can 
be  checked  by  hand  in  situations  where  the  linearised  elastostatic 
equations  can  he  solved  explicitly. 


4 


Observe  that  C(X)  S  Svm  (T^S,  T^B')  ,  the  positive  definite 

■  pos  X  X 

synunetric  linear  transformations  on  T^S  . 

Assume  we  are  given  a  smooth  stored  energy  function  W  defined  on 

oairs  (  AX')  where  C  S  svm  (T„B,  T„  S)  .  For  ^  S  C  .  the  stored  energy 

■  pos  X  X  - 

of  3  is  I  .WCX,C(X]  )dVfXl ,  wliere  C  is  the  Cauchv-Oreen  tensor  of  3 

and  dV  is  the  volume  element  in  B  .  The  fact  that  V.'  depends  on  3 

only  through  the  ooint  values  of  C  is  usually  called  material  frame 

indifference.  (See  Truesdell  and  Noll  [ 1965j  and  Marsden  and  Hughes 

[1978]  for  discussions.)  Since  C  is  a  function  of  F,  we  shall, 

T 

by  abuse  of  notation  also  write  W(X,F)  for  W(X,F  F) . 

The  first  Piola-Kirchhoff  stress  tensor  P(X,F)  is  defined 

9 

by  P(X,F)  =  ;rcW(X,F),  the  partial  derivative  of  W  with  respect 
0  r 

to  F.  Thus,  P(X,F)  e  L(T^6J1'^)  .  The  second  Piola-Kirchhoff  stress 

tensor  S(X,C)  is  defined  by  S(X,C)  =  -Ir^VCX.C),  so 
— — (j  c 

* 

S(X,C)  £  sym(T^B,  T^B)  .  From  the  chain  rule,  one  has  the  relationship 
P(X,F)-G  =  4s(X,C)-  ' fV,  *  G^F] 
for  all  G  e  L(T^8dR'^)  .  ' 

For  finite  dimensional  vector  spaces  V,W,  the  bilinear  map 

L(IV,V)  X  L(V,W)  -*•  IR;  (A,B)  -^race(A  B)  defines  an  isomorphism 
* 

L(V,1V)  =  L(iV,V).  If  V  =  W  is  an  inner  product  space,  this 

* 

isomorphism  identifies  syTn(V,V)  with  s>ti(A,\'1.  Using  this  identifi¬ 
cation,  we  get  PrX,F)  €  L(1R^,  T^B)  ,  S(X,C)  £  sym(T^B,T^81  and 
PfX,F)  =  S(X,C)‘’F^,  or  P  =  SF^  for  short, 

"TTie  mass  density  does  not  appear  in  our  formulas  as  we  arc  building 
it  into  the  definitions  and  work,  for  example,  with  hodv  force 
per  unit  volume  rather  than  per  unit  mass. 


o 


Let  A(X,F)  =  denote  the  elasticity  tensor.  Identifying 

L(LlT^S JR'^)  ,  L(T.^B,  IR"^)*)  in  the  usual  way  for  second  derivatives 
with  B(L(T^B,  IR"^)),  the  space  of  bilinear  maps  of 
L(T^B,  IR'’)  L{T^B,  IR"^)  to  IR  ,  A(X,F)  determines  a  symmetric 
bilinear  map  on  L(T^8,  IR'^^ . 

The  second  elasticity  tensor  CfXjC)  is  similarly  defined  to 


be  —  =  "* 
®  3C  - 


3“W 


evaluated  at  (X,C)  ,  so  may  be  regarded  as  a 
svTTunetric  bilinear  map  on  sym(T^5,  T^B) .  The  chain  rule  gives 


2A(X,F)  =  CCX.CI'CF^H  *  H^F.  F^G  *  G^F)  +  S(X,C)*(H^G  +  G^H) 


The  following  two  assumptions  will  be  made  in  the  first  two 
papers  of  this  series. 

(HI)  The  undeformed  state  is  stress  free;  i.e.  P(X,I)  =  0, 
or  equivalently,  S(X,I)  =  0,  where  I  is  the  identity. 

(H2)  Strong  ellipticity  holds:  there  is  an  e  >  0  such  that 

A(X,I)*(v®5,  V  ^  >  ell^i'r llvll" 

for  all  C  t  T.^^B  and  v  eiR"^,  where  v  0  ^  G  L(T  8,  R"^)  is  defined  !>\ 

'  '  (v  0  (V)-vf  l\  ,1 . 

The  classical  elasticity  tensor  is  defined  by  C(X)  =  2C(X,I) , 
so  c(X)  is  a  symmetric  bilinear  mapping  on  symuT^BjT^B)  to  R;  at 
({)  =  I,  we  identify  T^B  and  R"’  since  x  and  X  coincide.  By 


A(X,I)-(G,H)  =  ;|^IX)-(G  +  G^,  H  + 


(HI)  , 


6 


I 

! 

;  Regarding  ACX.I)  L(L(T^R,T^8') ,  L(T^8,  and 

C(X)  e  L(svt7i(T^B,  T^8).  sytn(T^S,  T^8))  ,  this  reads 

:A(X.I)-G  =  C(X)*(G  +  G^)  , 

or,  if  G  is  symmetric, 

A(X,I)-G  =  c(X)(G)  . 

By  (H2) ,  solvability  of  the  lineariced  equations  of  elastostatics  can 
be  determined  by  the  Fredholm  alternative  (see,  e.g.,  Marsden  and 
Hughes  [1978]). 

We  shall  let  B;  8  IR'^  denote  a  given  body  force  (.per  unit 
volume)  and  t:  38  -  a  given  surface  traction  (per  unit  area). 
These  are  dead  loads;  in  other  words,  the  equilibrium  equations 
for  ({)  that  we  are  studying  are: 


(E) 


DIV  P(X,  F(X))  +  B(X)  =0  for  X6  6 

P(X,F(X))-N(X)  =  T(X)  for  X  S  38 


where  N(X)  is  the  outward  unit  normal  to  38  at  X  €  36  and  DIV  P 
is  the  divergence*  of  P(X,F(X))  with  respect  to  X. 

Let  L  denote  the  space  of  all  pairs  fB,T)  =5,  of  loads  (of 
class  IV^  “’P  on  5  and  38  )  such  that 


f  B(X)dV(X}  +  T(X)dA(X)  =  0 
■'8  -*38 

‘Recall  from  above  that  P(X,F(X))  €  L(]R-^,  T^6) 
then  P(X,FfX))v  defines  a  vector  field  on  8  ; 
defines  the  vector  field  DIV  P  by  (DiV  P)-v  = 


For  any  v  fc  , 
its  divergence 
DIV  (Pv1  . 


i.e.  the  total  force  on  B  vanishes,  where  dV  and  dA  are  the 


respective  volume  and  area  elements  on  6  and  3B.  Using  the 
divergence  theorem ,  observe  that  if  the  pair  is  such  that  (E) 

holds  for  some  5  S  C,  then  (B,t)  S  L  . 

Throughout  the  paper,  the  group  SOf3)  =  {Q  £  LQR^ j  Q^Q  =  I 
and  det  Q  =  +1}  of  proper  orthogonal  transformations  will  play  a 
key  role. 

By  (HI)  ,  1)  =  Ig  (the  identity  map  on  S)  solves  (E)  with 
B  =  T  =  0.  By  material  frame  indifference,  4)  =  Q|8  is  also  a 
solution  for  any  Q  €  SO(3).  The  map  Q  h-  q|8  embeds  SO(3)  into 
C  and  we  shall  identify  its  image  with  SO(3) .  Thus,  the  "trivial" 
solutions  of  (E)  are  elements  of  SO(3) . 

Our  basic  problem  is  as  follows: 

(PI)  Describe  the  set  of  all  solutions  of  (E)  near  the  trivial 
solutions  S0(3)  for  various  loads  H  6  L  near  zero. 

Here  "describe"  includes  the  following  objectives: 

(a)  counting  the  solutions 

(b)  determining  the  stability  of  the  solutions 

(c)  showing  that  the  results  are  insensitive  to  small 
perturbations  of  the  stored  energy  function  and  the  load. 


This  section  is  devoted  to  the  geometry  of  the  load  space  L. 
Many  of  the  results  of  this  section  are  available  in  the  literature, 
but  we  gather  them  here  for  convenience. 

Before  beginning,  we  shall  recall  a  few  notations  and  facts 
about  the  rotation  group  S0(3j .  Let 

M,  =  LCR'^,  =  linear  transformations  of  to 
sym  =  {A  e  =  A} 

skew  =  {A  e  M, |A^  =  -A} 

We  identify  skew  with  so(3),  the  Lie  algebra  of  S0(3);  skew  and 
E"^  are  isomorphic  by  the  mapping  v  CE"^  v*  v  G  skew  ,  where 
v(w)  =  w  X  V.  If  V  =  (p,q,r)  relative  to  the  standard  basis,  then 


The  Lie  bracket  is  [v,w]  =  vX’W-wSv  =  -(v  x  w)  where 
V  K  w  €  M,  is  given  by  (v  ®  w) (u)  =  v<w,u>  .  The  inner  product 
is  <v,w)  =  4tTace  (v  w) ,  the  Killing  form  on  so(3).  Finally, 
exp(v)  is  the  rotation  about  the  vector  v  in  the  positive  sense, 
through  the  angle  Ivll. 

Now  we  turn  to  a  study  of  L  .  For  €0  and  IG  L  ,  we  say 


that  1  is  equilibrated  relative  to  ?  if  the  total  torque  in  the 
configuration  ■(>  vanishes: 
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I  4>(X)  X  BfX)dV(X)  +  [  (})fXJ  X  T{X)dA(X)  =  0 

where  I  =  From  symmetry  of  the  stress  one  sees  that  if 

I  =  (B,t)  e  L  satisfies  (E]  for  some  ))  €  C,  then  Z  is  equilibrated 
relative  to  ,J).  (An  easy  proof  uses  the  Piola  transform;  cf.  Marsden 
and  Hughes  [1978]). 

Let  denote  the  loads  that  are  equilibrated  relative  to  the 

identity  configuration  Ig  . 

Define  the  astatic  load  map  k:  L  x  C  -  M,  by 

=  B(X)  »4>(X)dV(X)  [  T(X)  ».D(X)dA(X) 

and  write  k(Jl)  =  k(il,  Ig) . 

We  have  actions  of  SO (5)  on  L  and  C  given  by: 


Action  of  S0(3)  0£  L:  Q)l(X)  «  (QB(X) ,  Qt(X)) 
Action  of  S0(3)  on  C :  Q<t>  =  Q°<t> 


Note  that  Q9.  means  "the  load  arrows  are  rotated,  keeping 
the  body  fixed."  We  shall  write  and  0^,^  for  the  S0(3)  orbits 

of  5-  and  <1.  Thus,  0  denotes  the  trivial  solutions  corresponding 
to  =  0. 


The  following  is  a  list  of  basic  observations  about  the  astatic 


load,  each  of  which  is  readily  verified: 
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(Ai)  I  is  equilibrated  relative  to  J  if  and  only  if 

k(£,^)  €  sym.  In  particular,  I  G  if  and  only  if 
k(2,]  G  sym. 

(A2)  (equivariance)  .  For  !.  G  L ,  D  t  C.and  .Q.,  € 

k(Q^i.Q2-M  =  QikC^.D)Q2^ 

In  particular,  k(,Q)l)  =  Qk(Jl) 

(A3)  (infinitesimal  equivariance) .  For  IGL,  6  G  C  ,  and 
G  skew, 

k(W^il,^)  =  W^k(jl,0),  k(£.W,(^)  =  -k(i,.5)  W, 

In  particular,  k(WJl)  =  Wk(2). 

Later  on,  we  shall  be  concerned  with  how  the  orbit  of  i.  G  L 
meets  L  .  The  most  basic  result  in  this  direction  is  the  following. 

5.1.  DaSilv.a's  Theorem.  Let  I  G  L.  Then  0„  L  /  0. 

Proof.  By  the  polar  decomposition,  we  can  write  k(/:) 

for  some  Q  G  SO (5)  and  A  G  sym.  By  fAZ)  ,  kfQi)  =  QkiU 

so  by  (Al)  ,  q£  G  L  .  ■ 
e 

Similarly,  any  load  can  be  equilibrated  relative  to  any  chosen 
configuration  by  a  suitable  rotation. 

The  concept  of  an  axis  of  equilibrium  deals  with  the  case  in 
which  o’  meets  in  a  degenerate  way. 


T 

Q  A 


=  A  G  .svm. 
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3.2.  Definition.  Let  i  E  L  and  v  E  TR'^ ,  llvll  =  1.  We  say  that 
-  e 

V  is  an  axis  of  equilibrium  for  Z  when  exp(6v!i  e  for  all 

real  9,  i.e.  when  rotations  of  Z  about  the  axis  v  do  not 
destroy  equilibration  relative  to  the  identity. 

A  number  of  useful  ways  of  reformulating  the  condition  that  v 
be  an  axis  of  equilibrium  are  as  follows. 

3.3.  Proposition.  Let  i  E  and  A  =  k(;l)  E  sym.  The  following 
conditions  are  equivalent: 

1 .  I  has  an  axis  of  equilibrium  v 

2.  there  is  a  v  S  IR'^  ,  llvll  =  1  such  that  vl  E 

3.  W  H-  AW  +  WA  fails  to  be  an  isomorphism  of  skew  to  itself 

4.  trace  A  is  an  eigenvalue  of  A. 

Proof.  1  =»  2.  This  follows  by  differentiating  e.xpiiv)'.  9 

at  9  =  0.  2  =»  1.  By  f-AZ)  , 

k(exp(9v)?J  =  'I  +  (9v)  +  ^'(9v)"  +  ...’kCi)  . 

Since  k(v’]  =  vkC’l  is  symmetric,  this  is  symmetric,  term 
by  term. 

2  =*  5.  Since  k(vi)  =  vA  is  symmetric,  vA  Av  =  0,  so 
W  K  AW  *  WA  is  not  an  isomorphism. 

3^2.  There  exists  a  v  SIR''’,  llvll  =  1,  such  that  vA  +  Av  =  0, 
so  kfvi]  =  vA  is  symmetric, 

3  4 .  Define  L  E  M.  bv  L  =  (trace  AjI  -  \.  Then  one  has 

.> 


the  relationship 
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(Lv)"'  =  Av  +  ^A 

(In  fact,  if  [u,v,w]  denotes  the  triple  product,  the  relationship 

[Bu,Bv,Bw]  =  det  B[u,v,w]  gives  [Au.v.w]  +  [u,Av,  v;]  +  [u.v.Aw] 

^  T 

=  (trace  A)  [u,v,w]  .  This  yields  (Lv)  =  i'A  +  A  v  ,  which 

gives  the  claimed  results  for  symmetric  A).  Therefore, 

Av  +  vA  =  0  if  and  only  if  Lv  =  0  ,  i.e.,  v  is  an  eigenvector 
of  A  with  eigenvalue  trace  A  .  ® 

3.4.  Corollary.  Let  l  £  and  A  =  k()l)  G  sym.  Let  the  eigen¬ 
values  of  .A  be  denoted  a,  b,  c.  Then  i  has  no  axis  of  equili- 
brium  if  and  only  if  (a  +  b] (a  +  c) (b  +  c)  /  0. 

Proof.  This  condition  is  equivalent  to  saying  that  trace  A 
is  not  an  eigenvalue  of  A.  B 

3.5.  Definition.  We  shall  say  that  'I  G  is  a  type  0  load  if  I 
has  no  axis  of  equilibrium  and  if  the  eigenvalues  of  A  =  k(£) 
are  distinct. 

The  following  shows  how  the  orbits  of  type  0  loads  meet  L^. 

3.6.  Proposition .  Let  IG  L^  be  a  type  0  load.  Then  0^ 
consists  of  four  type  0  loads. 

Proof.  We  first  prove  that  the  Sn(,A)-orbit  of  A  in  M_  under 
the  action  Q-QA  meets  sym  in  four  points,  "elativc  to  its  basis 
of  eigenvectors,  we  can  write  A  =  diag(a,b,c).  Then  C  ^  '  '^yn 
contains  the  four  points 


diag(a,b,c) 

diag(-a,-b,cj 

diag(-a,b,-c) 

diag(a,-b,-c) 


(Q  =  I) 

(Q  =  diag(-l.-l,l)) 
(Q  =  diag(-l,l-l)) 
(Q  =  diag(l.-l,-l)) 


These  are  distinct  points  since  (a  +  b) (a  +  c) (b  +  c)  ^0.  Now 
suppose  a,  b  and  c  are  distinct.  Suppose  QA  =  S  S  sym.  Then 

■y 

S“  =  A~ .  Let  be  an  eigenvalue  of  S  with  eigenvector  u.  . 

2  9  ">  T  ■> 

Then  S  u^^  =  yTu^  =  A“u.  ,  so  pT  is  an  eigenvalue  of  A“ .  Thus, 

2 

as  the  eigenvectors  of  A“  with  a  given  eigenvalue  are  unique, 
u^  is  an  eigenvector  of  A  and  is  the  corresponding  eigenvalue. 

Since  det  Q  =  +1,  det  S  =  det  A,  so  we  must  have  one  of  the  four 
cases  above. 

By  equivariance ,  sym  =  ^  =5ym  is  a  set  consisting  of 

four  points.  Now  0^  >  ^o  it  suffices  to 

show  that  k  is  one-to-one  on  0^^.  This  is  a  consequence  of  the 
following  and  (A21 . 


3.7.  Lemma.  Suppose  A  €  sym  and  dim  ker  A  ^  1.  Then  A  has  no 
isotropy :  i . e.  QA  =  A  implies  Q  =  I . 


Proof.  Every  Q  ^  I  acts  on  IR'  by  rotation  through  an  angle 
9  about  a  unique  a.xis,  that  i.s,  about  a  line  through  the 
origin  in  .  Now  QA  =  A  means  that  Q  is  the  identity  on  the 
range  of  A.  Therefore  if  0^1  and  QA  =  A,  the  range  of  A 
must  be  0  or  1  dimensional,  i.e.  dim  ker  A  J.  ■ 
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Finally  in  this  section,  we  study  the  range  and  kernel  of 
k;  L  M. 

j 

3.8.  Proposition.  1.  ker  k  consists  of  those  loads  in  for 

which  every  axis  is  an  axis  of  equilibrium. 

2.  k:  L  -*•  M.  is  surjective. 

j  - - - 

Proof.  1.  Let  i  G  ker  k.  For  W  €  skew,  k(W)L)  =  WkCi)  =  0 

so  S  by  3.5  every  axis  is  an  axis  of  equilibrium.  Conversely, 

if  Wii,  €  L  for  all  W  G  skew,  then  k(WS,)  =  Wk(.2,j  is  svmmetric 
e 

for  all  W;  i.e.  k(?,)W  +■  Wk(i)  =  0  for  all  W.  From 

(Lv)^  =  Av  *  vA,  where  A  =  kill),  and  L  =  (trace  .A)I  -  A,  we 

see  that  L  =  0.  This  implies  trace  A  =  0  and  hence  A  =  0. 

To  prove  2  introduce  the  S0(3)  invariant  inner  product  on  L: 

<Z,l>  =  <  B(X)  ,B(X))dV(X)  +  I  <  T(X)  ,T(X)  >dAl,X)  . 

^8  •'38 

T 

Relative  to  this  and  the  inner  product  trace  (A  Bl/2  on  M,,  tlio 
T 

adjoint  k  :  M,  L  of  k  is  given  by 

k^lAj  =  (B,t)  where  BfXj  =  AX-G,  tfX)  =  AX-G. 

and 

G  =  h  AX  dVfX)  *  f  AX  d.AfXl'^/l'dV  '  dA  . 

^'8  ^  i5  -  -'S  '38 

If  k^(A)  =  (0,01  then  it  is  clear  that  A  =  i).  It  follows  that 
k  is  surject i ve.  B 


I 
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3.9.  Corollary.  1.  ker  k  is  the  largest  subspace  of  that  is 

S0(3)  invariant . 

2.  k|  (ker  k)'*'  ;  (ker  k)'*'  M,  is  an  isomorphism. 

Let  j  =  (k|  (ker  and  write 

Skew  =  j  (skew)  Sym  =  j  Csy™) 


These  are  linear  subspaces  of  L  of  dimension  3  and  6  respectively. 
Thus  we  have  the  decomposition; 


S0(3)  invariant  pieces- 


L  =  Skew  ®  Sym  ®  ker  k 


L 

e 


corresponding  to  the  decomposition  M.  =  skew  •©  sym; 

1  T  1  T 

U  =  4(U  -  U  )  +  j(U  +  U  ) . 


.Vote:  Skew  and  L  need  not  be  orthogonal. 
-  e 


i4.  Equivalent  Reformulations  of  the  Problem 


Define  C  L  by  Pfo)  =  (-DIV  P,  P-\')  i.e. 


P(P)(X)  =  C-DIV  P(X,F(X}],  P(X,FCX)) -.veX)) 


so  the  equilibrium  equations  (E)  become  =  I  .  From  material 
frame  indifference  we  have  equi variance  of  b;  i.e.  J(Qc)  =  Qb(i), 
Standard  Sobolev  estimates  show  that  i  is  a  smooth  mapping  (see, 
for  example,  Palais  [1968]).  The  derivative  of  t>  is  given  by 


D$((j))‘U  =  (-DIV(A*Vu) ,  (A*7u) ‘NJ 


and  at  =  !„  this  becomes 


D$(Ig)«u  =  (-DIV(c*e) ,(c-e)-N) 


1  T 

where  e  =  ^^iVu  +  (Vu)  [. 

If  D<J>(Ip)  :  T-  C  L  were  an  isomorphism  we  could  solve 

“  8 

$(0)  =  i  uniquely  for  i  near  Ig  and  I  small.  The  essence 

^  •  u  n.?,rT  ^  •  •  u-  vsince  P(S0(3])  =  0 

of  our  problem  is  that  D<P(Ig)  is  an  isomorphism; 

Define  C  =  ("u  S  T.  C'  u(0)  =  0  and  7u(0)  €  svm}  .  From 
sym  y  . 

o 

(H2)  and  linear  elasticity,  we  have; 


4.1.  Lemma.  D<I>(I„)IC  :  C  L  is  an  isomorphism. 

-  ■  B  '  svTTi  svm  e  - * - 
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The  connection  between  the  astatic  load  map  k :  L  NL  and  J 
is  seen  from  the  following  computation  of  k-'?. 

4  2 .  Lemma.  Let  v  €  C  and  let  P  be  the  first  Piola-Kirchhoff 
stress  tensor  of  i  .  Then 

kfSC:))  =  '  pdV 
'3 

This  follows  by  an  application  of  Gauss'  theorem  to 

kCirbll  =  I  f-DIV  P)  X»  X  dVCXl  +  '  fP  X'i  SX  d.Ai  X'l  . 

■  8  '*38 

This  should  be  compared  with  the  astatic  load  relative  to  the 
configuration  J  rather  than  Ig;  one  gets 

k(0(b)  =  I  S  dV 

'B 

which  is  symmetric,  while  k('?(i))  =  k f T (b  1  , 1,,)  need  not  be. 

o 

To  study  solutions  of  $(<{))  =  I  for  1)  near  the  trivial 

solutions  and  Z  near  a  given  load  it  suffices  to  take 

S  Lg.  This  follows  from  Da  Silva's  theorem  and  equivariance  of 

Let  be  regarded  as  an  affine  subspace  of  C  centered  at 

I„.  Let  ?  be  the  restriction  of  T  to  C  .  From  the  implicit 
8  sym  ^ 


function  theorem  we  get: 
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Later  we  shall  show  how  to  compute  D“F(0)  in  terms  of  D$(Ig) 
and  C. 

Now  we  are  ready  to  reformulate  problem  (PI). 

(P2)  For  a  given  S  near  cero,  study  how  0^  meets  the  graph 

of  F  for  various  I  near  . 

Problems  (PI)  and  (P2)  are  related  as  follows.  Let  o  solve  (L) 
with  ,?.  €  L  and  Q  be  such  that  5  =  Q(J)  S  Then  $(5)  =  Q2.  , 

so  the  orbit  of  Z  meets  the  graph  of  F  at  ^(i)).  Conversely,  if 
the  orbit  of  Z  meets  at  $(o)  =  QH  ,  then  (+i  =  Q  solves  (E) . 
We  claim  that  near  the  trivial  solutions,  the  numbers  of  solutions 
to  each  problem  also  correspond.  This  follows  from  the  next  lemma. 


4.4.  Lemma.  There  is  a  neighborhood  U  of  I g  ij).  ^svm  that 

0  G  U  and  Q(|)  S  U  implies  Q  =  Ig. 


Proof.  Note  that  C  is  transverse  to  C\  at  Ip  and  Ip 
-  syra  I g  o  o 

has  trivial  isotropy.  Thus,  as  S0(3)  is  compact,  0  is  closed,  so 

S 

there  is  a  neighborhood  of  Uq  of  Ig  in  such  that 

Q| 8  G  Up  implies  Q  =  I.  The  same  thing  is  true  of  orbits  passing 


through  a  small  neighborhood  of 


I  by  openness  of  transversality  and 


compactness  of  S0(3) .  ■ 


If  0^  meets  in  k  points  =  ^(0^),  i  =  l,...,k  then 

•i.  are  distinct  as  4>  is  1-1  on  a  neighborhood  of  I,,  in  C 
1  o  sym 

If  this  neighborhood  is  also  contained  in  U  of  4.1,  then  the  points 


:o 


Qi  are  also  distinct  by  -i.i.  Hence  the  problems  IPI'I  and  iPdj 

are  equivalent. 

In  connection  with  the  action  (Q,A)  h-  i^A  of  S0(31  on  M,  we 
shall  require  some  more  notation.  Let 

Skew  (A)  =  ^  skew 

and 

Sym(.A)  =  •4(-'^  ^  sym 

be  the  skew  symmetric  and  symmetric  parts  of  A,  respectively. 

We  shall,  by  abuse  of  notation,  suppress  j  and  identify  Sym 
with  sym  and  Skew  with  skew.  Thus  we  will  write  a  load  1  €  L  as 
i  =  (,A,n)  where  A  =  k(?,)  €  M.  and  n  €  ker  k;  hence  2  €  precisely 
when  A  €  sym.  The  action  of  SO  (5)  on  L  is  given  by 

Qe  =  (QA,Qn). 

Using  this  notation  we  can  reformulate  problem  (Pd)  as  follows: 

(P3)  For  a  given  =  (A^^.n^^j)  €  near  cero,  and  I  =  (A,n) 
near  Find  Q  €  SO (31  such  that 

SkewlQ/'  F(Sym(QA)  .Qnl  =  0. 

Define  the  rescaled  map  F :  IR  <  -*■  Skew  by 


(3. da) 


(3.db) 


:i 


F(\,2)  = 

Since  F(0)  =  0  and  DFCO")  =0,  F  is  smooth.  Moreover,  if 
F(Z)  =  +  ...  is  the  Taylor  expansion  of  F  about 

zero,  then  F(\,:i)  =  ...  . 

In  problem  (F.)  let  us  measure  the  sice  of  Z  by  the  parameter 
i  .  Thus,  replace  =  <.  for  Z  near  zero  by  I)(b)  =  '^Z 

for  A  near  cero.  This  scaling  enables  us  to  conveniently  distinguish 
the  size  of  .i  from  its  'orientation'.  In  the  literature  I  has 
always  been  fixed  and  \  taken  small.  Here  we  allow  I  to  vary  as 
well.  Thus  we  arrive  at  the  final  formulation  of  the  problem. 

(P43  For  a  given  =  (A^.n^)  •=  L^,  for  I  near  and  \  small , 

find  Q  e  S0(5)  such  that 

Skew(QA)  -  'F(\,  SymiQA) ,Qn)  =  0. 

The  left  hand  side  of  this  equation  will  be  denoted  H(\,A,n;Q) 


or  H(,\,Q)  if  A,  n  are  fixed. 


We  shall  begin  the  analysis  by  giving  an  (almost  trivial)  proof 


of  one  of  the  basic  theorems  of  Stoppelli  ’  1938) 

5.1.  Theorem.  Suppose  I  G  has  no  axis  of  equilibrium.  Then  for 
A  sufficiently  small,  there  is  a  unique  i  S  and  a  unique  Q 

in  a  neighborhood  of  the  identity  in  S0f3)  such  that  ^  =  Q*^i>  solves 
the  traction  problem 

$(<)))  =  Aj?,  . 

Proof.  Define  H;  IR  x  S0(3)  -*•  Skew  by 

H(X,Q)  =  Skew(QA)  -  AF (A ,Sym(QA) ,  Qn) ) 

where  I  =  (A,n)  €  =  Sym  ■©  ker  k  is  fixed.  N'ote  that 

D,H(0,I)-W  =  Skew  ( WAJ  =  *  AW).  By  Proposition  3.3,  this  is 

an  isomorphism.  Hence,  by  the  implicit  function  theorem,  H(,\,Q)  =  0 
can  be  uniquely  solved  for  Q  near  I  S  S0(3)  as  a  function  of  A 
near  0  S  IR.  ■ 

The  geometric  reason  "why"  this  proof  works  and  the  clue  to 
treating  other  cases  is  the  following. 


The  only  other  comt)  lete  proof  in  F.nglish  we  know  of  is  given  in 
Van  Buren  1.1968',  although  sketches  are  available  in  (irioli  1963  , 
Truesdell  and  Noll  '1965^'lnd  Wang  and  Tn;esdell  IP'S'  .  Dur  proof 
rather  different;  the  use  of  the  map  F  avoids  a  series  of 
licated  estimates  used  by  Stoppelli  and  Van  Buren. 


IS 

comp 


3  .  - 


Lemma.  A  load 


},  6  has  no  a.\is  of  equilibrium  precisely 

when  L  =  (d  T^O.,  .  In  particular,  if  I  has  no  axis  of  equili- 

brium,  then  C,  intersects  L  transverselv  at  I 
-  i  -  e  - 

Proof.  The  tangent  space  to  at  I  G  is  "  . Wil'lV  *5  skew; 

and  the  pro.iection  of  this  into  the  complement  Skew  to  is 

Wi  ’i(WA  AW)  where  A  =  k(i)  .  The  result  then  follows  from 
part  3  of  3.3.  ■ 


We  have  shown  that  there  is  only  one  solution  to  I'fp)  =  M.  near 
the  identity  if  ,1  is  small  and  has  no  axis  of  equilibrium.  How 
many  solutions  are  there  near  the  trivial  solutions  S0(5)?  As  we 
shall  see,  this  problem  has  a  non-trivial  answer  which  depends  on  the 
t>'pe  of  1.  We  analyte  the  simplest  case  here.  Recall  from 
definition  3.. 3  that  a  load  ?,  S  is  said  to  be  of  type  0  if  i  has 
no  axis  of  equilibrium  and  if  A  =  kdl  has  distinct  eigenvalues. 

Loads  with  no  axis  of  equilibrium  occur  amongst  other  types  of 
loads  classified  in  the  next  section,  and  StoppelH's  theorem  5.1 
applies  to  them.  However,  the  global  structure  of  the  solutions 
("global”  being  relative  to  S0(3))  is  different  for  the  different  types. 
For  type  0  the  situation  is  as  follows. 

5.5.  Theorem.  Let  S  be  of  type  0.  Then  for  '•  sufficient  !>' 
small,  fft)  =  „  has  exactly  four  solutions  t,  and  :  ,  in  a 

neighborhood  of  the  trivial  solutions  S0(3)  C  C  (see  Figure  5  . 


Proof.  By  3.6,  0,,  meets  in  four  Doints.  By  5.1,  in 

vx-Q  e 

a  neighborhood  of  0  in  L,  (?,  ,  meets  ^  in  e.xactly  four  points, 

\x.o 

the  images  of  5^,  f-,.  S-  and  '5^,  say.  Thus  problem  (P2)  has  four 
solutions.  By  the  equivalence  of  (PI)  and  (P2) ,  so  does  (PI).  ® 


SO(3) ;  the 
trivial  solutions 


Let  A  =  and  =  {Q|QA  €  sym}.  From  tlie  proof  of  i.b 

we  see  that  3^  is  a  four  element  subgroup  of  sols')  isomorphic  to 
Z-,  ®  Z^.  By  our  earlier  discussions,  the  elements  5^  are  obtained 
from  by  applying  rotations  close  to  elements  of  S^.  In 

particular,  as  '\  0,  the  solutions  converge  to  the  four 

element  set  (regarded  as  a  subset  of  C) . 

For  i  sufficiently  close  to  the  problem  ?(;)  =  X  .1  will 

also  have  four  solutions.  Indeed  by  the  openness  of  t ransversa 1 ity , 

0^^  will  also  meet  .V  in  four  points.  In  other  words,  the  picture 
for  tv-pe  0  in  Figure  2  is  stable  under  small  perturbations  of  1^. 

Ne.xt  we  study  the  stability  of  the  four  solutions  found  in 
Theorem  5.5.  This  will  be  done  under  the  h.vpothesis  that  the  classical 
elasticity  tensor  in  stable.  IVe  introduce  the  following  condition; 

(H5)  Assume  there  is  an  n  >  0  such  that  for  all  e  •=  SvtdIT^S,  , 

e(e}  =  4<^(X)(e,e)  >  nilell',  il'll  =  pointwise  norm 

(£(e)  is  the  stored  energy  function  for  linearized  elasticity!. 

Because  of  difficulties  with  potential  wells  and  dynamical 
stability  in  elasticity  (see  Knops  and  Wilkes  (19~I'  and  Ball,  Knops 
and  Marsden  (ID'S]  1  we  shall  adopt  the  following  "energy  criterion" 
definition  of  stability. 

3.4.  Definition .  A  solution  d  of  M.-i)  =  1.  will  be  called  stable 


if  i 


is  a  local  minimum  in  C  of  the  potential  function 


:b 


V ,  i  ?)  =  Wi  ijdv  - 
j 

5 

r  r 

where  B(Xj  •DiX')dV[X)  +  ■  t  f  X)  •  :  (X )  dAlX)  =  trace  k  (  ?. ,  t  ; 

>  5  ^ 

If  0  is  not  stable,  its  index  is  the  dimension  of  the  largest 
subspace  of  vectors  u  tangent  to  C  at  D  with  the  property 
that  ;  decreases  along  some  curve  tangent  to  u.  (Thus,  index  0 
corresponds  to  stability.) 


5.3.  Theorem .  Assume  iHIi-CHSl  and  let  i he  as  in  5.3.  For 
sufficiently  small,  amongst  the  four  solutions  given 

by  5.3,  exactly  one  is  stable;  the  others  have  indices  1,  3,  and  3, 
Suppose  i  is  a  solution  approaching  Q  €  ^  '\  -  0.  Then  t 

is  stable  if  and  only  if  QA  -  trace(0A)I  €  s>m  is  positive-definite. 
In  general,  the  index  of  5  is  the  number  of  negative  eigenvalues 
21  QA  -trace  (QA) 1  . 


Proof.  Let  ~  C  solve  M:)  =  ■  ■  ,r  ■  •  Then  is  a  critical 

point  of  V,  ,  .  Consider  the  orbit  c\  =  'Qv,-.  Q  -  Sdi.iV  of  ; 

''■n 

Its  tangent  snace  decomposes  T,  C  as  follows; 

•0 


T.  C  =  T.  t\  -f-  IT.  0. 

•  ()  '0  *0  'D  ‘II 


first  consider  V,  restricted  to  IT.  0.  .  I 

'  ‘U  ‘11 


at  :  in  the  direction  of  a  -  iT.  f.  ■'  i« 

-It  *o 


'(•cond  leri-ativ' 


.:-f  ■! 


•  u  ■ 


QIS 


this  becomes 


[  c(X).(e(X),e(X)jdV(X), 


where 


1  T 

e  =  — (7u  +  (7u)  ).  This  is  larger  than  a  positive  constant  times  the 

7 

L"  norm  of  e,  by  [H3^ .  However,  since  u  is  in 
2  2 

Heir,,  >  (.constant)  lluir.  by  Korn's  ineouality  (see  Fichera  '  1972  (). 

l~  H 

By  continuity,  we  have 


D-y, 


'^0 


(u,u)  ^  lluir 
H 


if  u  is  orthogonal  to  0.  at  and  \  is  small.  This  implies 

i>0  0 

that  tn  is  a  minimum  for  V, ,  in  directions  transverse  to  C 

0  Pq 

(.Actually  one  can  see  that  is  a  local  minimum  in  the  topology 

I 

of  C  on  (T  0  )"  by  using  Tromba's  '29"6'  version  of  the  .Morse 

’O  ^0  ' 

lemma. ) 

Next,  consider  V,,  restricted  to  C  .  Bv  material  frame 
0  '0 

indifference,  W  is  constant  on  c*,.  and  so  as  i  must  be  a 

critical  point  for  V.  restricted  to  0,  .  it  is  also  a  critical 

'*^'0  "0 

point  for  \Z„  =  I  restricted  to  C  (where  <ltl  =  It 

suffices  therefore  to  determine  the  index  of  '0.  at  The 

result  is  now  a  consequence  of  continuity  and  the  limiting  case 
■  -*  9  given  in  the  following  lemma  about  type  0  loads. 


5 . P .  Lemma .  Let  I  he  type  0  and  let  A  =  k ( f ) .  Then  S  ^ , 

regarded  as  a  subset  of  C  equals  the  set  of  critical  points  of 

C\  .  These  4  critical  points  arc  nondegenerate  wit.n  indices 

8 

9 ,  1  ,  2  and  5 ;  the  index  of  <)  is  the  number  of  negative  «■  i  genv.a  L.'e 


of  )A  - 1  r.ace '  (IM  I . 


Proof.  E-irst  note  that  L  =  IT,  SOljVr  since  3(51  I„1  nas 
■  e  1  ^  0 

o 

kernel  T  S0f5)  =  skew,  lias  range  L  and  is  se If-aJi oint .  Thus 
o 

Q£  S  /.  if  and  onlv  if  1  i  T  S0l5).  It  follows  that  QH  S  L  if 
e  T  ^  e 

j  0 

and  only  if  Q  is  a  critical  point  of  id  .  .Recall  that 
elements  of  3^  =  {Q  <=  St'('5):Q£  £  are  symmetric.  ) 

To  compute  the  index  of  ,,  d  at  Q  6  S  we  compute  the 

ig 

second  derivative 


t=0 


■l(iV-Q3 


N’ow 


£(WQ)  =  trace  kf  i.W'Q)  =  trace[kf ,Oi'a 


- 1  ■’ 

=  trace  iAQ  w~  ]  =  tr.ace  [WQ.A]  ' 


because  Q  =  Q  •  Tliis  ouadratic  form  on  skew  is  represented  bv  the 
elonent  O.A  -  trace(0A)I  s>it!  a.s  is  seen  fror.  v.X  +  .Xv  =  '  I.\’ i  with 
■A  replaced  by  OA  and  tr.ace  !  vw  i  '  dvw  .  ''siny  tlie  rc'^resentat  ions 
for  fQA-  given  Proposition  j.n,  namcly 

Jiagla.b.ci,  di ag < -a , -b  ,  c 1  .  di ag ( -a ,h , -c 1  and  diagi a , -h  , -c 1 


one  checks  that  all  four  indices  occur.  P 


Remark.  This  lemma  is  ,i  ipoci.il  case  of  the  ceneral  problem 


to  study  tb.o  critic.il  points  ■'!'  ::neir  '■'unct  ;  (in.;  1  s  ,ui  orbits  oi' 


a 


f 


■'ll 


representation  of  a  Lie  group.  This  situation  will  arise  again  in 
our  analysis  of  the  other  load  types,  of.  Frankel  [19651  ana 
Ramanu j  am  _ 1 9  6  9 1 . 


iO 


Class i fication  of  Orbits  in  M.. 

- o 

The  purpose  of  this  section  is  to  classify  orbits  in  under 

the  action  (Q,A)  v*  qa  of  S0(3J  on  M-  by  the  way  the  orbits 
meet  sym.  It  is  enough  to  consider  orbits  0^  of  elements  of  sym 
by  the  polar  decomposition.  We  begin  by  recalling  a  result  already 
proved. 

b.l.  Proposition  (,T\'pe  0).  Suppose  A  £  svcn  has  no  a.xis  of  equilibrium 
and  has  distinct  eigenvalues.  Then  0^  sym  consists  of  four  point. s 
at  each  of  which  the  intersection  is  transversal. 

This  was  Proposition  3.6.  (.Another  proof  of  this  is  given  below.} 

We  shall  let  the  eigenvalues  of  A€  sym  be  denoted  a,b,c. 

Using  the  terminology  from  53,  we  say  that  A  has  no  a.xis  of  eoui librium 
when  (a+b) (b+c) (a+c)  ^  0;  i.e.  a  +  b  +  c  ^  a,  b  or  c ,  and  in  this 
case  0^  intersects  sym  transversely  at  A. 

6.3.  Definition.  We  shall  say  A  is  of  type  1  if  A  has  no  a.xis 
of  equilibrium  and  if  exactly  two  of  a,b.c  are  equal  and  non-tero 
I  say  a  =  bi^c,  a^O). 

6.3.  Proposition .  I f  A  is  type  1,  then  sym  consists  of 

two  points  (each  with  no  a.xis  of  equilibrium)  and  an  IRT^  ^  '  eacn 

point  of  which  has  one  axis  of  equi  librium'i  . 

Before  proving  this,  we  give  a  number  of  lemmas  of  genera! 
utility.  If  fc  IRP"  is  a  line  through  the  origin  in  ,  let  gi 
be  the  rotation  through  angle  ~  about. 


» 


.  -I  .  Lemma  . 


,  is  an  embeddiniJ  of  HP”  onto  iOtS' 


Proof.  It  is  clear  that  I  Q,  is  a  one-to-one  map  of  PP" 
into  S0(3).  Since  QT  =  I.  Q.,  =  Q,  =  Q.,  so  it  maps  into 


SO  1.3]  sv-m. 


F.very  S0!3]\I  is  a  rotation  through  some  angle  -  aiout 


some  a.xis 


.  If  such  Q  also  is  symmetric  then  it  has  three 


indenendent  real  eigenvectors.  Hence  r 


0.5.  Coro  1  lary .  The  orbit  of  the  identity,  0,,  meets  s>m  i.a 
one  point  (I)  and  IRP”  =  (S0(3)  sym'A  [. 


6.6,  Lemma.  Let  A  €  sntti  with  dim  ker  A  i  1  and  suppose  that 
Q  S  SOfj]'  I  and  QA  €  sym.  Then  0  =  0,,  for  some  line  ■.  invariant 
under  A  ,  and  in  narticular  fP  ^  s)!!!. 


^roof.  We  can  suppose  Q  d  i.  There  is  a  unique  :up  to  .s  i  cn 

unit  vector  .\  G  IR'^  Such  that  Qx  =  x.  .since  QA  -3  .-^ym,  we  have 
7' 

=  AQ  ,  so  QAQ  =  A.  Thus  QAx  =  .\x .  so  \x  =  cx  for  a 

con.st.int  c.  Lach  of  Q  and  \  Ic.aves  \  =  x  ,  the  ortiiogonal 

complement  of  x.  invariant,  an.!  \  is  not  identicalLv  zero  on 

Lot  S  =  OA  t.  syn,  so  3“  =  .\”  .  Since  Q'\’  is  a  rot.it  ;on 

thi.s  implies  S'-\'  =  C_  A'V  giving  0=1  or  ^  =  Q,  where 

'  ’  !  X  1 

«.  ;x'  is  the  line  through  x  .  Then  m  sym  .as  in  Lemma  i' .  i  .  ■ 


-<1 


It  follows  t 
QA  S  3y;n,  then 
A  and  Q  can  be 


hat  if  dim  Ver  \  1  iin! 

QA  =  AQ,  so  as  0  is  both  orthogonal  and  symmetric, 
simultaneously  diagona li red. 


Proof  of  0.1.  If  .A  has  u*jtinct  eigenvalues,  its  eigenvectors 
are  unique,  so  0  is  either  the  identity  or  is  a  rotation  by  "  about 
one  of  the  eigenvectors.  ■ 

Proof  of  “'.5.  Suppose  0  =:  a  =  b  f  c ,  and  let  w  he  an  eigenvector 
correspo.’iding  to  the  eigenvalue  c.  Let  V  be  the  plane  orthogonal 
to  .»  ;;  :<  t.te  eigenspace  with  eigenvalue  a.  As  Q,  A  can  be 

-  i.tul'' I'v.'ou;  ■  .iiagonalited  and  Q  is  a  rotation  by  m  (excluding 
'  =  :  .sv  ,.i’.e  vither  Q  =  or  Q  =  Q,  for  i  a  line  in  V. 

f  :'rmer  case,  Q.  ,^.|A  has  eigenvalues  (-a,-a,c)  so  has 
•O'  >'f  ouui  i  ibriiut.  In  the  latter  case ,  Q.,.A  has  eigenvalues 

.1,-1, -Cl,  i.'i  w  1 -s  an  axis  of  equilibrium,  {see  3..S).  B 

b.~.  Corollary.  The  in  Proposition  6.5  is  a  right  coset  of 

the  subgroup  S^,  of  all  rotations  about  w;  in  fact  =  s|(^. 

where  is  any  fixed  line  in  V,  the  plane  orthogonal  to  w. 

Proof.  MP^  =  {Qj I V}  and  we  have  the  easily  verified 


identity 


=  exn( -iv  ..j, 

' t/:  ■  "0 

where  c,  rakes  an  an?le  -  with  in  the  •Kisiti'/e  sense  in  . 

.  -  ' 

These  lemmas  also  enable  us  to  handle  the  next  t>'pe . 

6..S.  Definition .  IVe  shall  say  A  is  of  type  i  ;f  \  has  no 

axis  of  equilibrium  and  all  three  of  a.b.c  are  equal  iand  so  i  . 

6.9.  Proposition.  I f  A  is  type  d,  then  J',  '  s)!!!  consists  of 
one  point  (A)  and  an  IRP“, 

Proof.  This  is  immediate  from  6.5.  9 

.Motice  that  each  point  of  the  W  has  a  whole  circle  of  axes 
of  enuilibria;  namely  Q^A  has  as  axes  of  equilibria  all  vectors 
orthogonal  to  i.  The  eigenvalues  of  Q,,A  are  a, -a, -a. 

Types  0,1,  and  2  exhaust  all  symmetric  matrices  with  no  axis 
of  equilibrium,  and  it  is  easy  to  check  from  the  above  that  any 
symmetric  A  with  dim  ker  A  c  I  lies  on  the  SO (5) -orbit  of  a 

type  0,1,  or  2  matrix.  From  now  on  we  shall  call  these  orbits,  or 

any  representat i ves  of  them,  tv-pe  0,1,  or  2 . 

Finally  we  turn  to  the  remaining  A's  with  an  axis  of  equilibrium 
that  is  not  already  on  an  orbit  of  t>'pe  0,1,  or  2. 

We  say  A  is  type  5  if  dim  ker  A  =  2  ai'.d 
f  A  =  0. 

.  I  f  A  j_s_  type  5,  then  0^  ■Aym  cvin'':sts 


6.10.  Definition. 


Proof.  S  =  0.4  S  svTTi  implies  S~  =  A“  anJ  so  again  S  =  ;^\, 
as  in  0.6  even  though  possibly  .\'V  =  0  .  In  this  case  Q  could 
be  any  rotation  about  ’iC.xi.B 

,-Vll  the  foregoing  information  can  be  summarited  as  follows: 

6.12.  Theorem.  The  50(3)  orbits  in  .M _  fall  into  five  distinct 
types  according  to  the  way  in  which  they  .meet  sv-m  (  sec  Table  1  below.)  . 
Furthermore,  for  A  S  svan,  =  vQ;QA  -3  sym;  consists  of  1  {Q, } 

for  all  I  invariant  under  A  (and  hence  ~  i>'m)  except  in  the 
cases  (1)  dim  ker  ,A  =  2 ;  th:’‘n  also  contains  the  rotations 

through  any  angle  about  the  eigen-axis  of  A  corresponding  to  the 
non -zero  eigenvalue; 

(2)  A  =  0;  t hen  =  S0{3)  .  See  Table  2  below. 


Remarks 

1.  Table  1  highlights  the  fact  that  having  an  axis  of  equilibrium 

or  not  is  not  an  invariant  of  the  S0(3)  action  on  L  .  This  means 
there  are  equilibrated  loads  having  an  axis  of  equilibrium,  but 
which, wiien  rotated  globally  by  a  certain  amount  to  another 
equilibrated  load, no  longer  liave  one. 

2.  Thus,  b^  Theorem  5.1,  we  get  existence  of  solutions  to  tl'.e  traction 
problem  for  all  t vpie s  of  astatic  loads  except  3,1.* 

3.  The  notion  of  tvpe  can  be  nulled  back  from  M_  to  L  with  a 

.5 

little  care,  as  we  sec  below. 


*In  particular  the  occurrence  of  0  solutions  !\v  Stoppelli  in  ‘yne  1 
is  seen  to  be  due  to  a  neglect  of  the  full  rot.ation  group  i  see 
Section  S) .  Our  results  are  also  consistent  witli  those  of  Ball 
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b . 13.  Definition . 

By  analogy  with  our  definition 

=  (Q  e  SO (5)  I  QA  e  s>Tn}  ,  A  e  M. 

which  we  applied  when  A  is  of  type  0  ,  let  u.s  now  write 

S,j  =  {Q  S  S0(3)  ;  0^  €  L^'  ,  i  ^  L  . 

From  equivariance  of  k  we  clearly  have; 

6.14.  Lemma.  . 

■Vote  that  the  map  -*•  0^0  sym:  Oi-^  QA  is  an  embedding  for 
types  0,1,2  but  not  all  t>pes  3,4,  because  of  the  isotropy. 
Pulling  back  to  L  ,  we  see  that  Q»—  0?,  is  an  embedding  S,  - 
if  kCHj  is  of  type  0,1  or  2,  so  we  can  refer  to  ?  as  being  of 
tvpe  0,1,  or  2  according  as  k(i)  is.  On  the  other  hand,  if  k(?.1 
is  of  tvpe  3  then 

either  (a)  0  L  = 

or  (b)  Oj  3  =  two  disjoint  circles  in  *  ker  k 

F'inally,  if  k(l)  is  of  type  4  then  0,  3  ker  k  ^  and  any 
SO (31  orbit  C,  is  allowable. 


ure  l•lgenvaluus 


[  Mil  s  \  r>.i  1  i  nt  f  r  sec  t  1  vJii  nt  tlic  iu*bit  ot  A  with  syin. 


types.  These  loads  are  all  pure  traction,  with  B  =  0. 


Type  L .  Rotacion  by  180® 
about  one  of  the  horizontal 
axes  produces  an  equili’-- 
rated  load  vith  no  axis  of 
aquilibrium. 


Type  2 .  Any  horizontal  axis 
is  an  axis  of  equilibrium; 
vertical  axis  is  not  an  axis 
of  equilibrium.  Rotation  by 
130®  about  the  vertical  axis 
gives  an  equilibrated  load 
with  no  axis  of  equilibrium. 


Type  3  (o.)  .  The  load 
itself  admits  a  circle 
group  of  summe tries 
about  the  axis  i-- 
which  is  thus  an  axis 
of  equilibrium. 
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Type  2  (  o) .  The  load  is 
not  syyanetric,  but  the 
astatic  load  remains 
constant  under  rotation 
about  Che  axis  — which 

is  thus  an  axis  of 
equilibrium. 


Type  ■i  .  The  astatic 
load  is  zero:  all  axes 
are  axes  of  equilibrium. 


.oad  types 


'•D 


3  T .  The  Bifurcation  Equation  and  its  character 

According  to  the  formulation  ihl'i  of  our  problem,  .<c  wi>h  to 
solve  the  equation  H(\,A,n;Q)  =  0  for  Q,  where 

Hl'.\,A,n;Q)  =  Skewi'.lA)  -  '-F  (  \  ,Syn(QAl  ,  An  )  , 


lA,n)  is  near  iA,,,n,J  i  and  is  small.  Define  the  ri  ght  -  i  nvar ;  ant 

i  U' 

vector  field  X.  on  S0(3)  by 

X  (Q.)  =  Ske-.v(QA  ].Q: 

i.e.  right  translation  of  SkewCPA^)  €  so(3)  =  1^30(3)  to  T|^S0(3). 
Likewise,  we  shall  regard  H  as  a  right- invariant  vector  field  ci 
30(3)  depending  on  the  parameters  ,\,A,n  by  setting 

X(\,A.n;Q)  =  HC- .A.n  ;Q) -Q 


Thus , 


X|0,A,,,n,,:Q)  =  iQ). 

0 


Final  Iv,  note  t!i.-it 


he  zero  set  of  \ 


=  ■  0  SOI  31  Sivew'‘./-\  ,  = 

I' 


What  s 


for  vari''us  t  vpes  was 


1  S 


r^o  P-  i 
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.  1 .  Lemma. 


Suppose 


sym  is  of  any  type.  Then  *^or  0 


T  S 
Q  A. 


i.WQjW  S  skew  and  ''‘QA^  + 


Proof.  The  second  equalitv  is  clear  for  anv  ,  because  (Qj : 

■  ■  d 

'VQr—  skew  1,'VQAyj'Q  .  For  the  first  one,  t;ie  inclusion  C  iirjt.ediately 
follows  by  differentiation  of  X.  (QJ  =  0  in  0.  Equality  then 

■ '5 

follows  by  a  dimension  count;  recall  from  3..S  that  v  t-  v  gives 

an  isomorphism  from  the  space  of  axes  of  equilibrium  for  A  (not 

necessarily  of  unit  length)  to  the  IV  €  skew  such  that  KA  +  A'.V  =  i).  ■ 

T 

Recall  that  W  WQA^^  ♦  corresponds  to  the  linear  trans¬ 
formation  t  race  (QA.  J I -OA  ,  under  tlv'  i  somo'r'h  i  of  •-'ke'..  =  SOfXl 

'0  ii 

with  When  Q  ^  S  ,  QA  is  symmetric,  so  tliis  t  rans  format  i  on 

N) 

is  symmetric  relative  to  the  Killing  form  on  This  remark 

and  ".1  proves  the  next  lemma. 


Lemma . 


:)u:^nose 


::r  inv 


ne . 


Then  at  each  point 


ot  5^  ,  the  ran^e  of  DX,  iQl:  T^SOlSi  T^SO'5')  is  the  orthoconal 

—  A,,  -  0  Q  - ^ - 

complement  ot  T ^  . 


Next  we  recall  a  general  context  tor  the  bifurcation  of  vector 
fields  that  will  be  applied  to  our  situation  icf.  Reeken  'IA'3'V  Let 
M  and  be  manifolds  and  A;  M  <  A  —  TM  a  smooth  vector  field  on 
M  depending  on  the  parameters  i  €  A.  We  seek  the  teros  of  X.  For 
''  =  '  ,  suppose  the  zero  set  S  of  X  is  a  known  smooth  compact 

submanifold  of  M.  Assume  that  M  carries  a  Riemannian  metric  and 
that  for  X  €  S,  the  range  of  D^^X(x,\q)  is  the  orthogonal  complement 
of  T^S.  The  normal  bundle  E  of  S  trivializes  a  neighborhood  U 
of  S.  For  each  x  €  U,  let  P  ;  T  M  T  S  ,  .  be  the  orthogonal 
projection  to  the  fiber  over  t(x)  ,  where  t;  E  ^  S  is  the 

projection.  By  the  inverse  function  theorem,  there  is  a  unique 
section  5,  :  S  -  E  such  that  P^^X( 0,  ( x)  , \)  =  0  for  x  €  S  and 
i.n  a  neighborhood  of  fassume,  for  example,  that  *  is  comf^act  ). 

Let  Xfx,!]  be  the  orthogonal  projection  of  X(x,')  onto  the  tangent 
space  to  the  graph  of  t.  at  a  point  x  on  the  graph.  Thus,  X  C  x , ' 
is  a  vector  field  on  the  graph  of  f.  and  finding  its  zeros  is  clcarl}’ 
equivalent  i  for  small  '■  ]  to  finding  zeros  of  X.  We  call  the  e  paation 
X(x,'‘l  =  n  on  the  graph  of  i-  the  bifurcation  equation.  Since  s 
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and  the  graph  of  are  di ffeomorphic  under  r,  ,  we  can  equally 

well  regard  X  as  a  vector  field  on  S.  This  reduction  of  the  nroblen 

is  often  known  as  the  Lianunov-Schmidt  method. 

The  above  procedure  may  be  applied  to  our  vector  field 

X(A,A,n;Ql  with  parameters  (\,.A,n)  and  variable  .k  =  £  50(51  = 

Assume  \  is  near  tero  and  fA.n)  is  near  a  load  where  A, 

is  of  arbitrary  t>’pe .  Thus,  there  is  a  unique  section  of  the 

normal  bundle  to  S.  determined  by  the  Liapunov-Schmidt  procedure 
■^0 

as  described  above.  Let  TfA.A.n)  denote  the  graph  of  \  ^ 

and  let  X(A,A,n;Q)  be  the  orthogonal  projection  of  X  to  the 
tangent  space  of  7  at  Q.  Thus,  X  is  a  vector  field  on  7. 

As  above,  we  may  also  regard  X  as  a  vector  field  on  5,  . 

■^0 

The  rest  of  this  section  is  devoted  to  proving  that  the  essential 
part  of  X  is  a  gradient.  In  the  general  context  above,  if  X  is 
a  gradient,  then  so  is  X  since  the  orthogonal  projection  of  a  gradient 
vector  field  to  a  submanifold  is  the  gradient  of  the  restriction.  This 
simple  version  does  not  directly  apply  to  our  situation  as  X  need 
not  be  a  gradient  vector  field  on  S0(3).  However,  the  "second  order" 
Taylor  approximation  X,  of  X  will  be. 

To  state  our  gradient  results,  recall  that  in  ,'4  wo  defined  the 
quadratic  function  G:  -*  skew  to  be  the  second  order  term  in  the 

Tavlor  o.xnansion  of  F  about  0.  Thus  F --Gf.'l  ... 

h 

where  G  is  a  quadratic  function  of  1.  The  .appropriate  second  order 
approximation  to  the  vector  field  X  will  thus  be  defined  !yv 


^s  k  ew  ( ,  I  (4 . 


X,i  '■  ,A,n: 'j1 


skew ( OA 1 
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Let  X,  be  the  second  order  approximation  of  the  vector  field  X  on 

S  obtained  by  the  Liapunov-Schmidt  orocedure.  Tlius ,  X^fQ)  is  the 

'^0 


orthogonal  projection  of 


.3.  Theorem.  Suooose 


gradient  vector  field  on 


f[Q)  -  < 


X.,  onto  the  tangent  space  -or  Q  €  s  . 

■■  '0  0 


is  of  arbitrary  type.  Then  X.,  is  a 
S  .  In  fact,  X,  =  -grad  f,  where 


<7UQ,c(7uy)dV 


J8 

and  =  Dofl  i.e.,  is  the  uniaue  solution  in  C 

-  Q  ^  ^  0^  ’  - ’  Q  - - -  s\Tn 

of  the  Lineariced  equations  with  load  . 


Recall  that  the  pairing  between  loads  ■.  =  ■'B,')  and  configurations 
(or  displacements)  is  given  by 

^  B(X)*C(X)dV  *  '  'X-;  \  dV  =  trace  k '  C , ' 

:.6 


and  physically  represents  a  potential  for  ti'.e  working  of  the  loads. 

_  T  .  T  . 

Observe  that  if  I  ^  L  ,  then  ^  ''.O'i.v  =  trace  '  AO  i  ^  tracolAO  )  =  {  i  ,QI 

e  '  0 

for  all  Q  €  S0(,3)  . 

Remark .  In  the  second  term  of  X^  and  f  we  can  replace  ■',! 
bv  !.  .  However,  the  difference  is  higher  order,  so  ,  is  sufficient 
for  subsequent  applications. 


Proof  of  ~ .3.  We  shall  show  that  X,  is  a  gradient  field  on 
SOI  31  which,  by  the  remarks  follwoing  “.3.,  is  sufficient. 

Wc  proceed  in  two  parts.  l.et  us  first  stiow  that  X^'Oi  ;s  the 
gradient  of  (  on  all  of  .S(tl3i  . 


7.4 

Lemma . 

Let 

-  L 

and 

A  =  ki 

fij.  . 

Define  a  vector  field 

on 

S0(3) 

X.lQI 

—  Skew ( QA) 

•Q  as 

above 

and  the  man 

of 

SO  (3) 

IR 

■  iQj 

=  (  i 

,Q^I^  , 

.  Then 

\'^  =  f. 

Proof.  Two  si.-nple,  but  useful  observations  are: 
If  B,  W  fc  M, ,  with  W  S  skew,  then 


<  B,IV>  =  (  skew  B,1V  ^ 


and 


If  B  G  M. ,  I  G  i  and  t  G  C  >  then 


<  i,B*>  =  <  > 


To  prove  ".4,  we  compute  as  follows: 


d^(Q)-(WQ)  =  <  (WQ)’^Ig  ) 

=  <  (KOj^,  k(y„))  by  (; 
=  <  (IVO)^.A  ) 

=  (  K  ^ , OA 


-  W ,  skew  iOX  i,'  bv  '  I 
=  WO,  in.\yn  ' 

=  -'WO,  \^o 


T 


4" 


This  deals  with  the  first  term  of  X^.  To  deal  with  the  second 
term,  we  need  a  special  case  of  Betti's  reciprocity  theorem: 

7.5.  Lemma.  =  <  (WO)?,^,u^)  for  and  i'WQ)  2,^  S  S\ti  . 

This  is  a  trivial  consequence  of  svmmetrv  of  DOfl^j  i.e.  of  the 

o 

elasticity  tensor.  It  is  also  proved  in  standard  references;  for 
e.xample,  see  Truesdell  and  Noll  ’  1965';  p.  525. 


To  prove  2’.  5.,  we  shall  al.so  need  to  calculate  the  .■second  derivative  of  the 

skew  component  of  •>;  i.e.  of  F(P)  =  Skew:  k(^ (o))j  .  Surprisingly, 

this  second  derivative  depends  only  on  the  classical  elasticity 

tensor  d  ■  Recall  from  52  that  we  regard  c  as  a  linear  map  of 

1  T 

sym  to  itself  and  that  we  write  e  =  ^rCTu+C^ti)  }. 

7.6.  Lemma.  Let  f:  C  skew  be  defined  by  =  Skew  . 

Then  F(Ig)  =  0,  DF(Ig)  =  0  and 


D“F(I  )(u,ul  =  2  Skew  (I  7u-cfe)dV  =  -2  Skew  k(?.  .u) 
t-  Jg  u 

where  "  -niV(c(e))  and  =  c(e]*N  ,  Identifying 

skew  with  31'^,  this  becomes 


-D“FfIg)  lu,u) 


b  '<  u-  dV  +  !"  T  '  n  dA 
J  U  j  LI 

6  38 


IS 


Proof.  By  Lemma  4.2,  K' ?}  =  SkewT*  Pd\' '  where  I’  is  the 

'  '  8 

first  Piola-Ki rchhof f  stress  tensor.  We  have 

5 

/*  ''P  " 

F(I„)  =  0.  Also,  DF(r„)-u  =  skew  .  dV  =  skew  e*e  d\  = 

6  S  )6jF 

tp  ° 

as  C’e  is  symmetric  and  since  — ( I  )  =  c.  To  comnute  ii  - , 

jt-  B 

shall  need  to  use  the  fact  that  S  is  symmetric.  Write  I'  =  : iiui 

use  the  product  rule  to  obtain  DpPlF)*7u  =  Tu-SiPl  +  Fn^SiF''.‘u. 

Thus,  as  ^(r„)  =  0, 

0 


\ 

D-P  (Ig)  •  (7u,7v]  =  7u*DpS^lg)  •7v  *  7vDpS(Fi -70  ^  D"  SI  I  g)  '  '  ~u  ,  7v\ 


T 

.Now  D-S(Ip)*7u  =  D^S(rc,)*(7u  +  7u  )  =  o*e  and  Dustin')  is  symmetric, 
r  o  ^  0  r  o 

SO 


D-FfIJ*(u,u)  =  skew  !  d:T(I„)  (7u,':u)dV 
o  j  F  b 

6 


2  skew  (1  Tu'Cte)  JV 


8 


Finally,  this  equals 


-2  skew 


U.  dV 


'55  u.  d.'V  } 
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by  the  divergence  theorem,  so  the  last  statement  follows.  T 


CxampK*.  For  i  hono^eneous  isotronic  T.aterial, 


:  I  e  ’  =  '•  I  trace  e)  I  -  2ue 


where  e  =  [Tu  *  l''u)  j/2  and  u  are  the  Lame  moduli. 


D“F(I_)iu,uj  =  2  Skew  ('  -'\7u*rtrace  (.ul]I  2u7u-e;d\' 
®  'S 


=  2  Skew  [trace(7u)  ]7u  +  u  .u'7u;dV'  M 

's 


Let  us  ne.xt  see  what  ".b  says  about  the  quadratic  term  G  in 

the  Taylor  e.xpansion  of  F.  For  i  £  C  we  have  the  identitv 
^  ^  svm 


Fl.f)  = 


where  L  is  the  orthogonal  projection.  Thus,  as  DF  and 

DF  are  zero  at  ami  0  resnect  i  vel >•,  and  P  , 

o  '  c  5  5 


D“F(U)(u,v)  =  D-FCOUD-fdJ-u.D.fd^.l-vl 
b  b  b 


Hence  for  i  S  and  writing  u,  =  1  ,  wo  have  the  identitv 

e  -  S' 


Gill  =  D“F(di(u,,ud 
b  <• 


50 


-G(;i)  =  2  Skew  [  [  b  ®  u  dV  +  [  t  ®  u  dAj 

"  j  *  J  '' 
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=  2  Skew  kUb.T).u^) 


where  b  =  -DIV(c*  (e)  ,) ,  and  x  =  c(e,)-N  and  e,^  =  [.u,  )  ]  '2  . 

However,  these  last  equations  say  exactly  that  (b.i)  =  I,  and  so  we  get 


-  =  Skew  k(,l,u^) 


Completion  of  the  proof  of  ".5.  The  derivative  of 
i  T 

Q  H-  <  Zq,  ^  Up)  in  the  direction  WQ  is  given  by 


<  a Q.,  T(woj^Up)  *  (  > 


<  ,  C'''Q)  >  (by  Betti  reciprocity,  ".S'l 


.(  w  .k(Q;,,^,Up)> 


by  (2) 


•  <  )V  ,  Skew  k  ,Up,)  >  by  (1) 


-<  WO  ,  Skew  k(Q;|^,Up)Q  ) 


<  WO,  4^1(QC^).0  ) 


b  V  I  5 )  .  I 


28 .  Bifurcation  Analysis  for  Type  I 

We  now  discuss  the  solutions  of  the  basic  equation 

H(:\,A,.M;0)  =  Skew(OA)  -  \F(A,S>Tn(OAl  ,Qn)  =  0  i  1' 

for  the  load  i  =  (A,n)  near  a  load  =  (A^.n^^J  having  an  axis  of 

equilibrium  and  of  type  1,  and  for  ^  near  d  .  We  shall  also  obtain 

the  stability  of  the  solutions  and  finally  we  shall  compare  our  results 

with  those  of  Stonpelli  [1958].  For  type  1  we  need  to  do  a  bifurcation 

analysis  on  the  circle  S  corresponding  to  the  degenerate  zero  set 

'■^0 

of  H  when  =  0  and  I  =  .  The  analysis  has  some  features  in 

common  with  the  papers  of  Hale  [1977]  and  Hale  and  Taboas  [J9S1], 

Without  loss  of  generality  we  can  assume  that  A^^  =  dia?  ,  ■  a , -c  1 
where  0  a"  /  c‘.  Thus,  from  §6  the  set  3.  of  zeros  of  Skew  fOA]  for 

■'o 

n  S  30(3}  is  given  explicitly  by  the  following  two  points  and  circle: 


S  =  '"diagC  1 . -1  ,-l)  diagf -1 , 1 1)  /  U  C 

■■'ri  0 


Cl 


where 


,  -y  »  . 

Q  ■=  I  y  x  U  1  I X  =  cos  ,  y  =  sin 
(1  1  / 


The  loads  corresponding  to  the  two  points  are  A^.^  =  diag  ia,a,ci  and 

+•  ★ 

A^^  =  diag  I  -  a  ,  -  a  ,  c  1  . 

From  ".3,  we  are  led  to  study  the  critical  points  i^f 
T  ■  T 

fnyi  =  (  .  ,1}  ^  '(1’^^  '*A ''  itnm  the  divergence 


theorem 


hat 


Thus  the 


'  o 

where  =  D?(n  ^  and  e^  =  [-u^  *  fTu^j“]/:  .  Thus  the 

function  f  is  computable  from  linearited  elasticity  alone,  which 

leads  to  the  curious  observation  that  our  "second  order"  nonlinear 

elasticity  here  involves  no  more  data  than  linear  elasticity,  but 

merely  processes  the  information  in  a  different  way.  Writing  Q= 

X  -y  0  \ 

y  X  0  j  as  in  (2),  f  becomes  a  noivnomial  of  degree  2  in 
1  0  II  1  / 


(x.y).  Write  the  two  terms  of  f  as 


f(,Q)  =  fix.y)  =  ;h  r  b^x  w  b,y)  +  *  ^3-'’  *  '‘r'  *  ""s-'  * 


f  (x,y)  =  T-f(.x,y) 


14) 


which  defines  the  numerical  constants  by,b^,b^  and  a^ . a^.  Sext 

define  new  parameters  by  writing 


and  letting 


(3' 


f  fx,y)  =  XjX-  *  a,xy  >  u.y-  *  UjX  *  a.y  *  a,. 


Note  that  depend  on  our  parameters 

the  elastic  moduli  of  the  material.  Thus, 


as  well  as  on 


^2  =  =*2’  ‘5  =  ^*3 


‘4  "  ^1  *  ‘3  "  ^2  *  -‘t-  ‘i.  "  ‘  '‘i 


Replacing  ■  by  ,  '.vhere  Q  is  as  in  '2',,  effects  a 

rotation  of  the  x-y  plane.  Thus,  by  rotation  of  i  if  necessary, 
we  can  assume  =  i)  . 

Let  us  fix  a,  ,  1,  and  consider  the  bifurcations  of  zeros  of 


df*  ,, 

=  2( a.- a,  XV  + 
dr  0  1  • 

cl 

-  a , V  on  S  1 . e . , 

,  critical  points  of 

4"  ’ 

1 

on  with  a^ 

and 

a_  as  parameters. 

a 

Set  M  =  l(a  ,■'1  SR'  <  (a  ,a.f1  =  ,  the  manifold 

4  o  dr  4  :> 

of  critical  points  of  f * .  Indeed,  M  is  a  manifoll  and  can  be 
parametrized  bv  .t:  R  <  -*•  M  ,  pfu.i'l  =  (-dia.+u)  cos  -2('a,*a)S 

1  O  ,71 

">  1  2 

Denote  by  r:  R  x  S  -*■  R  the  projection  onto  the  first  tactor, 

3.1.  Lemma.  Set  L  =  [2(0, -a,)'  -  a  "-'t.'J'  -  10S,a  'a.' ( a, -a.i”  .  f; 

-  -  Id  aa  4slo 

If  a,  -  a.  ^  0  ,  then  r,  M  R“  is  a  proner  stable  nan  in  ^  a ,  ,  a.  1  -  space . 

—  1  j  -  - - - - - , -  j  j  • — ' 

and  its  set  of  critical  values  is  the  astroid  defined  by  =  () 

I  see  Figure  below)  . 


df‘ 
d  ■ 


Since  the  number  of  points  in  ”  '(ti  (i.c.,  the  zeros  of 

at  a  =  (a,,a_)l  is  a  constant  over  ■  i)  or  ■  0  ,  we 
4  a 


obtain 


S.2.  forollarv. 


df* 

d-' 


has  4  zeros  if 


■  ;) 


and.  has  2  zero; 


0 


Proof  of  Lemma  3.1. 


.1.J 


The  critical  set 
•:  iL.,9)  s  R  X  3^  j 

critical  values  of  tt 


I  of  t-jq:  K  <  -►  ]R“ 

sin-Q  +  a,  cos~0  +  a  =  0}  . 

.5 

can  be  parametriced  by 


is 

Thus,  the  set  of 


a ,  =  -2  (a,  -a,)  cos'^  9 
4  1  j 

=  2(a^-a,)  sin^  9 


fS) 


Since  2  consists  of  4  cusp  points,  4  fold  lines  and 
T’cll  is  1-1  ,  by  a  result  of  \'?hitney  (see  Mather  [1969]  or 
'k)lubitsky  and  -uillemin  [1973]),  one  knows  that  Ttop  is  a  stable  map. 

Eliminating  6  produces  the  bifurcation  set 

[2(4  -a  )]-^^  =  a  ‘  (9) 

1  J  4  o 


For  otj^  -  ct,  /  0  ,  (S)  describes  the  astroid  shown  in  Figure 


Figure  4- 


Next ,  observe  that  for  real  numbers  A,  R  and  C, 


•A  +  B  +  C  =  0  if  and  only  if  A'^  +  R'^  <■  =  3ABC 


by  virtue  of  the  identity  A’’  *  R'^  ♦  C'’  -  3ARC  =  'A  +  R  *  fi'A" 

AB  -  BC  -  CA).  Applying  (101  t<i  fOl  shows  that  .  ,  i .apiivalen 

d  f  d  ,  J/3  J/3 

~t ,  *  i.  -  d  I  a,  -  r,  )  =  -  3i£  ,  i,  (  d  (  I ,  -  i_  !  1 

A  .I  1  •)  A  .)  1  o 

gives  the  stated  conclusion.  B 


('uhinr  |■>l)th 


The  family 


of  functions  on  S  with  parameters 


i  »  - 

4  o 


jf* 

dO 


enjoys  a  universal  property.  Consider  a  perturbed  family 

df*  -  m 

-75—  +  g(.\,c,f)  ,  with  g(0,0,'5''  =  0  tor  e'>t;  -K  <  r<  .  To 

dy 

,  df* 

each  (''.,c),  denote  by  M,  ^  '(^3—  +  g  j  {  \ ,  c  ,  ,  i.  ,y 

the  "manifold"  of  teros. 


S.3.  Lemma ■  For  (\,c)  sufficiently  small,  the  sets  M._  _  are 

manifolds  and  there  e.xist  two  smooth  families  of  di ffeomornhisms 

ti,  :  1R~  IR”  and  'f,  :  M,  -*  '1  defined  for  , c  sufficientlv 

•  ,  c  -  \  ,c  ■  ,  c  -  - ^ 

small,  such  that  ''’T.  =  ■>,  ,  and  T,,  =  identity,  =  identity. 

- -  A ,  c  ■  X ,  c  -  0,0  -0,0 


1  1 

Proof:  For  '',c  sufficiently  small,  the  map  t.  IR  ■  S  -  M,  _  . 

^  ^  *  ,  c 

If'  ’  '■''‘3*'’'^  ^  ^  ’  ' 

defines  a  paranetrizat ion  of  M.^  ^  .  By  Lemma  S.l,  "’C.  11  ^  -  IR 

is  an  unfolding  of  the  proner  stable  map  .  Tlius ,  one  can  find 

- !  2 

di ffeomorphisms  i*.  ^  ,  f,  on  IR  >  S‘  .ind  F“  resprct  i  ve  1  >• 

such  that  i'- ('■^'-0,  )  =  .  Th ;  <  lemma  follc.ws  lettine 

A,c  ' ,  c  '  ,  c 

F,  =  F*.  -.:r^  .  ■ 

'  '  .  c  '  ,  C  '  ,  c 


how,  we  are  re.ady  to  ^t.ate  our  n-ain  result  on  the  number 
solution-'  near  type  I  loads.  Let  ■  =  •' ( c  depend  ^mnot.'ily  -ui  a 
parameter  ■:  in  Jl”'  ,  with  •  - 1  =  .  Recall  rliat  A 

Je  fined  by  ,  '  -  ,  a ,  ,  a,  and  .i_  bv  1  i  . 

i  -  >  ‘ 


3.4.  Theorem.  Let  oe  a  type  1  load  with  =  >a,-a,-c), 

0  ^  a~  c~  ,  =  0  ,  and  a,  ii  a,  .  Then,  there  exists  a 

_  -  1  o  - 

(smooth]  function  TiXiCl,  T(a,0)  =  T(a^,a_')  +  OC- J  defined  tor 
I,  ,  c  1  sufficiently  small  and  \  >  0  such  that,  the  traction 
problem  has: 

( i  1  four  solutions  for  the  load  1  ( c )  i_f  J. ( '■  , ')  0 

(two  of  them  near  )  ■ 

'■0 

(ii)  six  solutions  for  the  load  'c(c)  i_f  3.(4, c)  >  0 
(two  of  them  near  C  ]  . 

Proof:  The  bifurcation  of  zeros  of  .\  fcf.  on  C.  (=S^) 

'o 

is  the  same  as  finding  zeros  of  —  '*  —  gf'  .c,") 

=  ^+-rg  or  ^  +  g,  where  g(0,0,f]  =  0  .  Let 

ut".  -  “  ut*  ■“  ■“  •  ,  C 

be  the  family  of  diffeomornhisms  found  in  Lemma  S.3.  Take 
b  (L(c)]  b,{',(c]) 

AC'.c]  =  L’T.  ,  (  — a  ,  -  *  a.]  which 

‘  ,  C  '  -4  '  O 

has  the  desired  property.  1 

Next,  we  want  to  determine  the  "generic"  structure  of  the  bifurcation 

set  K  =  =  0:  in  (>,c)  space,  '•  >  0. 

Lor  m  =  0  ,  (oe,,a„)  =  .',(a,,a.)  ^  0  .  it  is  clear  that  K.  -  nH  . 

4  o  4  s 

Indeed,  our  traction  problem  has  two  so i ut ions  near  C  it 

'o 

■(3  ,a.l  (1  ,  and  four  solutions  near  C,  if  (.■i,,a_)  0  . 

1  a 

b|A(C;  b,\<~l 

Fbnr  m  =  1  ,  consider  k;  c  >— *  i — -  *■  3^,  — -  <■  a.l 

where  ilc)  is  the  line.ar  p.irt  of  .ic'  .  T!i  i represents  a  sTr:u'cht 
line  winch  we  assume  to  intersi-ct  tiie  astroid  t  ran  sx 'tso  I  y  if  ‘iiey 


I 


meet.  \'otice  that  <-  =  C'',  .<’ 


.5  tae  inverse  image  o! 


the  astroid  (defined  bv  equation  ('91),  under  the  man 
b^,cic)  '  b,Ucl 

hy  :  — y,  L— -i-  '  •  as  in  the 

\  _  C  4  D  ^  ,  C 


proof  of  Theorem  3.4).  Write  \c  =  c  .  Recall  that  i(cj  =  Afc)  +  0(  c") 


and  consider  the  map 


C  H-h,  f.i.c)  =  y,  - 

^  .  V  C  — 


0(a1. 


b.,.Afc) 


OC'.l)  . 


Since  the  astroid  is  bounded  and  y,  is  close  to  the  identity, 

\  ,c 


there  e.xists  an  interval  such  that  K, 


^c  S  K.  }  Ci  f-M,Ml 


for  \  >  0  ,  c  sufficiently  small.  .Applying  the  isotopy  theorem 
for  transversal  maps  (see  e.g.,  Hirsch  [lO^b])  to  the  family  h- 
through  Iiq  =  k  ,  we  conclude  that  the  bifurcation  set  'X  consists 
of  0,  2,  or  4  curves  with  slopes  k  ^(astroid)  (see  Figures  .5,5) . 
Thus,  for  e.xample,  by  choosing  c  ^  0  sufficiently  small,  and 
letting  \  ^  0  (i.e.,  consider  the  load  \x.(c)),  one  can  pass  from 

a  parameter  region  where  there  are  two  solutions  near  the  circle 
(four  in  all]  to  one  where  there  are  four  near  the  circle  (six  in 
all)  and  back  again  to  the  two-solution  region  (see  Figure 
Such  a  situation  is  not  dealt  with  in  the  analysis  of  Stoppelli  [1938]. 


For  m  a  2 


let  us  suppose  that  the  affine  tiap; 


b  A(c)  b,A(cl 

^ -  *  a  ,  — -  *  a.)  is  onto,  wiiere  again  Al'c) 

is  linear  part  of  (c)  .  Without  loss  of  generality,  we  may  also 
assume  that  b^ACcl  =  anJ  b^A(c)  =  c.,  ,  where  c  =(c^,c-,,a')  . 


Notice  that  K,^  _  =  {  (c^  ,c, ,Cj  ,c,  ,  z)  S'*'-;  is  the  inverse  image 


of  the  astroid  under  the  map  h^  _ ;  (CpC^l'^  .  , 


h^  A(  i 


>  ^  -s  ♦  - 


b,A(.'0 


*  a.  =  a. 


,  j''  ^  +  a  ,  ^  +  a^  )  .  Set 


•  .  'C^  =  c_  .  and  consider  the  man 


As  betore, 


IS  bounded  unitorml 


for  \  >  0  ,  c  sufficiently  small.  Applying  the  isotopy  theorem 

for  transversal  maps  to  the  family  h,  ^  through  the  affine 

isomornhism  h,,  ,  we  conclude  that  the  bifurcation  set  is  a 

0, 0 

cy  1  inder- 1  ike  set  along  the  t-axis  with  '-sase  a  cone  over  the 


astroid  in  o^,c.,  space.  The  first  order  approximation  of  this 
cone  is  given  by  the  cone  over  the  astroid  in  the  plane  =  1, 
centered  at  r-2a,,  -2a-1  with  "site”  4'a,-a,  (see  Figure?}. 

4  a  1  j  ■ 


7 


Figure  7 


X'exr  we  discuss  the  stability  ot  the  solutions  correspondin';  to 


loads  near  load  of  t>T)e  1.  This  can  be  determined  by  combining 
our  stability  results  for  t^.nie  d  .Theorem  5.3;  together  with 
well-known  stability  results  for  the  cusp.  We  make  the  same 
assumotions  as  those  in  Theorem  3.4. 


3.5.  Theorem.  Let  ^  =  diag la , -a , -c  )  ,  =  diag(a,a,c.l  and 

*  * 

Aq  =  diag(-a, -a,c)  as  above.  The  indices  of  the  bifurcating 

solutions  are  bo.xed  in  Table  5.  IRecall  that  stable  solutions  have 

inde.x  =  Q.~)  In  each  case  the  circle  represents  C  defined  by 

0 


equation  (2)  . 


•Note  that  two  stable  solutions  bifurcate 
c  >  a!  .  In  all  other  cases  the  solutions 
unstable . 


off  the  circle  when 
near  the  circle  are 


L 


64 


which  is  a  constant  [independent  of  x,  y) .  In  this  situation, 

=  ;t  =  0  and  so  our  theorems  do  not  apply. 

3.5.  Example ■  Consider  the  same  traction  problem  as  above,  but 
with  a  homogeneous  nonisotropic  hyperelastic  material  whose 
linearized  elasticity  tensor  is  given  by  cfu)  =  ^  iiag  ‘2  ■ 


In  this  case, 

u  (X)  =  TqI  0 

0 

-a 

0\-| 

O'  X  , 

Q=f; 

-y  0^ 

X  0  ; 

^  U  Vo 

0 

-c  /- 

1 0 

0  1  ^ 

where  c  ^(F)  =  F  +  diag  F  .  Then 

/ 2ax  ay  0  \  /ax  ay  0  \  \ 

ay  '2ax  0  )  ,  [ay  -ax  0  ,  ) 

V  o'  0  -:c  •  \  0  0  -c  / 

-> 

=  4a"x“  +  2a"y''  +  2c“ 

1  ^  ,  ,  , 
Hence,  1  =  8a  “  ^  0  ,  and  our  traction  problem  tn  1.1'  nas 

six  .solutions  [four  near  )  ,  with  stability  determined  by 


table  .5. 


\e.\r  -.ve  sluiil  aiscass  :u)v*  to  obtain  tiio  results  of  ftoppelli 
1958  as  a  special  case  of  our  analysis.  SVe  refer  the  reader  to 
the  statements  of  Stoppelli's  results  in  Grioli  '19b2,  p.  58'.  l.a 
this  approach  one  focuses  attention  on  bifurcations  that  occur  on  the 
circle  by  e.xamininy  what  happens  at  a  particular  location  on  the 
circle  and  ■  .  V.’e  c.an  .issume  this  coint  is  il.'i)  i.e.  r  - 

u':rh  no  loss  of  general  ity. 

First  of  all,  if  t-,  +  a.  ?!  i)  then  (1,0)  is  not  a  critical 
point  of  f * ,  so  there  are  no  solutions  near  (1,0).  He  may  a.ssume 
then  that  +  a.  =  0  and  then  the  Taylor  e.xpansion  of  f*  about 
9=0  becomes 


.  iba, 

o  1 


(liighor  order 


of  Grioli  19<i2  . 


i 

Case  3.  If 

+  :i- 

j> 

-  -r^  =  ''  and  \  , 

t  then 

df  * 

W'  " 

6  .  2 
T-t  >  ’ 

+  Old’i  and  so  t 

iiere  are  i  ,  ■  r  _  .i 

i  :ir  i  x'iiis 

: fold  point ] . 

This  is  Theorem 

1  on  p .  ."  s  J  r  irioi: 

l'Jf'2 

i 

Case  3.  If  -1^ 

'^3 

-  T  =  -  a  ■  - 

,  ‘  1  , 
but  t.  -  ^  *  1 

1  .1  ^ 

! 

.  d  f  *  4 

TT  "  3''h 

-  .  or,*.. 

so  there  ,ire  I,  or  . 

) 

solutions  'cusp  point).  This  is  Theor'cm  H  on  or  ii-bli  '.'.'''.1 


Furthermore,  it  we  express  our  constants  x.  i=a.  '  ;n  r  c  r-.-s  ot  t;ie 
elasticity  tensor  c  and  solutions  or  the  lineariced  probu'm 
usinii  f3)  above,  we  find  the  same  conditions  for  t'nese  ' 

as  is  given  on  p.  S"  of  Grioli  Idhd.. 

Thus  wo  recover  the  results  of  Stopnelli  on  type  1  loads,  \s 
was  explained  in  the  introduction,  however,  this  analysis  is  on!;.' 
local  on  tlie  circle  and  does  not  give  the  full  stor;.'  of  tlie 
bifurcation  picture,  even  in  this  c:ise.  I'iie  comph  '.e  hii'ui'cati  -i 
analx'sis,  including  stability. is  summariced  by  our  iMgurv  iS.; 

Table  3. 
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